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ABSTRACT 

Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense 
gas varies across the circumnuclear regions in starburst galaxies. 

Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2 — 1) and CS(3 — 2) lines in the 
circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected formaldehyde (H 2 CO) 
and methanol (GH 3 OH) in both galaxies. We marginally detect the isotopologue C^"^S. 

Results. We calculate column densities under LTE conditions for CS and CH 3 OH. Using the detections accumulated 
here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer 
model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the 
galactic centres. 

Conclusions. Using the rotational diagram method for CH 3 OH, we obtain a lower limit temperature of 14 K. In addition 
to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the 
dense gas as traced by CS (and CH 3 OH). We also estimate the quantity of the dense gas. We find that, provided 
there are between 10® and 10® dense cores in our beam, for both target galaxies, emission of CS from warm (T = 
100 - 400 K), dense (n(H 2 ) = 10®“® cm“®) cores, possibly with a high cosmic ray ionisation rate = 100 1 ^ 0 ) best 
describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further 
from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient 
CS formation, which means that gas that is less dense by an order of magnitude is required to replicate observations 
in some cases. 

Key words, galaxies: starburst - galaxies: star formation 


1. Introduction 


The formation of stars, in particular the most massive ones, 
is one of the most important processes regulating the evolu¬ 
tion of galaxies; to trace this process, it is crucial to identify 
the chemical nature of the gas involved and to determine the 
physical conditions that are more prone to lead to star for¬ 
mation. Massive stars usually form in large concentrations 
of dense, warm gas. These conditions are able to survive dis- 
ruptive winds and radiatio n from nearby newly formed stars 
( jZinnecker fc Yorke|2007 1. In particular, starburst galaxies 
are powered by exceptionally high rates of massive star for¬ 
mation, possibly triggered by mergers between galaxies rich 
in interstellar matter. Such events must be relatively short¬ 
lived, persisting only until the interstellar gas reservoir is 
significantly depleted. 

Subject to a proper interpretation, observations of 
molecules in starbursts can be used for many purposes: 
tracing the leftover reservoir of matter from the star for¬ 
mation process; tracing the star formation process itself; 
and determining the galaxy energetics through the i nflu¬ 
ence of young stars on the surroun ding environments (Ke- 
pley et al. 20141 Sliwa et al. 20141. More specifically, the 
detection of molecular star-forming gas is one of the most 


direct ways to measure the star formation rate and activity 
in a galaxy, allowing us to significantly improve our under¬ 
standing of galaxy formation and evolution. To date, the 
wealth of molecular data for at least the nearest galaxies, 
shows a chemical diversity and complexity that cannot be 
explained by a one-component, static model, and indicates 
how relative abundances between molecules may be able to 
provide insights into the physical distribution of the molec¬ 
ular gas and the energetics of these galaxies. 


While some observed molecular transitions may arise 
from a UV-dominated gas, many can only be explained by 
the presence of large reservoirs of dense, warm gas (possibly 
pre-processed o n dust grains) that are a part of the star for - 
mation process (|Gao & Solomon||2004| |Onishi et al.||1998|) . 
In dense regions (n(H 2 >10'* cm“'®), the galactic interstellar 
radiation held is unimportant for the chemistry since it oc¬ 
curs in conditions of high visual extinction. Molecules that 
are particularly enhanced during the star formation pro¬ 
cess (including the warm up phase that follows the birth of 
protostars) are therefore good tracers of the formation of 
massive stars in galaxies. 


From a chemical point of view, sulphur-bearing species 
have been extensively shown to be particularly enhanced 
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during massive star formation in our own Galaxy, as well as 
in ne arby galaxies ( Mauersberger et al.||1989 Martin et al. 
20061. Among sulfur-bearing species, carbon monosulphide 
(CS) appears as one o f the best tracers of very dense gas 
(Bronfman et al. TTo^ . Observationally, CS has been iden- 
tihed in several studies of nearby galaxies (Mauersberger 
et al. I [19891 [Martin et aT1|2006| |Bayet et al. |2008a[ |A1- 
adro et al.|2011 1. Following theoretical studies (Bayet et al. 


2008b 


|201ip, 

bFTD 


multi-line observations of the CS molecule 
in nearby (L ) < 10 Mpc) extra galactic environments were 
carried out (Bayet et al.| 20091. These studies show that 
despite the lack of spatial resolution it is possible to point 
to a reservoir of gas of relatively high density, ^ 10^ cm“^, 
traced by the low J (up to 4) transitions of CS and an even 
higher (> 10® cm“®) density component traced by the high 
J (up to 7) transitions. So far surveys have focused on CS 
detections in the centre of galaxies. To gain insight into the 
distribution of this gas in the environment surrounding the 
nucleus, galaxy mapping in CS was required. 


2. Target galaxies 

In this work, the spatial distribution of CS in two nearby 
galaxies, NGC 6946 and NGC 3079, is studied. We detect 
CS J = (2 — 1) and J = (3 — 2) lines in the central region 
and various offset positions in each galaxy. NGC 3079 is 
a nearby Seyfert 2 galaxy with a high inclination angle of 


pernovae, a lso with nuclear starburst activity (Engelbracht 
et al. ]T996t . We observe CS lines in the central point of each 
galaxy as well as at offset points staggered at increments of 
10” (Figure [B. At a distance of «19 Mpc and wS.S Mpc , 
we observe 92 pc arcsec“^ and 27 pc arcsec“^ in NGC 3079 
and NGC 6946, respectively. As a consequence, for NGC 
3079 we map out to a maximum of 2.8 kpc from our cen¬ 
tral pointing, while in NGC 6946 we map a smaller region, 
810 pc. Further information on both galaxies can be found 
in Table [l] 


Veilleux et al. (19941 discovered the line emitting su- 
perbubble in NGC 3079 extending 13” to the east of the 
nucleus. Within this bubble, they find extremely violent 
gas motions ranging over 2000 km s“^, which is also seen 
to the opposite, western side of the nucleus. In addition 
to this, unusual excitation shows that shocks may be im¬ 
portant and could be a significant contributor to line emis- 
sion a cross the bubble. From CO observations, |Koda et al.| 
(20021 show that the central 2 kpc region of the galaxy 
contains a smoothly distributed « 2 x 10® M(d of molec- 
ular gas. Also using CO observations, Leroy et al. (20111 
observe that NGC 6946 is also highly rich in molecular gas, 
both in the central region we observe and out into the de¬ 
tails of the spiral arms. Starburst activity in t he nucleus of 
NGC 6946 is more moderate than NGC 3079 (Turner & Ho 
19831. NGC 6946 has a small 15” bulge, within a larger 63” 
bar around it. Dense gas has been studi ed in the central 
region of this gala xy using tracers HCN (iSchinnerer et ah 
20071 and HCO'*' (Levine et al.||2008l. Levine et al. (2008|| 
study the central kpc and hnd that the density of the gas 
is roughly constant across the area studied and find that 
a plausible distribution is of dense clumps bound by self¬ 
gravity, interspersed throughout the lower density molecu¬ 
lar gas. However, they also state that it could be possible 
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Fig. 1: Maps of NGC 6946 in CS (J = 2 - 1) and CS (J = 
3-2) [top two] and the same for NGC 3079 [bottom two]. 
Units are arc seconds. Av = 8 km s“^ 


that the dense gas is homogeneously dis tributed about the 
galactic centre. Schinnerer et al. (20071 observed in great 
detail the central 50 pc of the galaxy and find that the 
HCN(1 — 0) intensity distribution in the very centre of the 
galaxy is centred around two peaks. They find that these 
peaks correspond quite well to the molecular gas as traced 


Article number, page 2 of 14 

































































































G. Kelly et al.: Mapping CS in starburst galaxies: Disentangling and characterising dense gas 


by C0(2 — 1). They state that one possible interpretation 
is that these gas peaks are the result of two spiral arms 
connecting in the central 100 pc. Although there have been 
detections of CS i n the centre of many nearby galaxies (e.g. 
Zhang et al.S2014 1 , a mapping and full analysis of its distri- 


bution in those galaxies has not yet been undertaken. The 
location and the characteristics of this gas away from galac¬ 
tic centres is crucial to our understanding of the processes 
and conditions that a necessary for massive star formation 
to occur. 


3. Observations 

The observations were carried out in August 2010 using the 
IRAM-30m telescope to include the CS(2 —1) and CS(3 —2) 
lines (i/ = 97.980 GHz and 146.969 GHz, respectively). In 
order to capture other species (methanol, formaldehyde) 
and isotopologues of GS, each observation was not centred 
upon the GS line. Instead a line survey from 95.5 GHz to 
99 GHz and from 144 GHz to 147.5 GHz was performed 
in order to detect as many chemical tracers as possible of 
the dense, star-forming gas. At the observed frequencies of 
GS(2 - 1) and GS(3 - 2), IRAM 30m has a HPBW of « 
25" and 17" with main beam efficiencies, Heff) of 0.81 and 
0.74, respectively. These beam sizes correspond to between 
approximately 450 pc and 650 pc (NGG 6946), and 1.5 kpc 
and 2.3 kpc (NGG 3079).. During observing, approximately 
every 2 hours, the pointing, focus and calibration were com¬ 
pleted on planets and evolved stars. The pointing error was 
estimated to be < 3". The observations were carried out 
under varying weather conditions. Tgys was between 100 
K and 300 K with medium-good weather conditions {T 225 
<0.15). The observations were done by wobbling the sec¬ 
ondary mirror with a beam throw of 180" in azimuth. For 
the observed spectra, the antenna temperature (T^*) was 
converted to main beam temperature {T^h) using T^b = 
(Feff/Heff)Ta*, where is the forward efficiency of the 
telescope (values ranging from 0.93 to 0.95). The channel 
width spacing was 4 km s“^ and has been smoothed down 
to a 8 km s“^. The EMIR EO and El receivers were used 
in configurations to optimise observing time. We combined 
this with the Wilma backend. The data reduction was car¬ 
ried out using the GLASS program, part of GILDAS soft- 

war43 


4. Results and analysis 

4.1. Results 


Maps of GS spectra for both galaxies are shown in Figure [l] 
Broadly, the strongest emission is seen in the central posi¬ 
tion (see Tables [^and [^for coordinates) of both galaxies, 
with two exceptions. In^NGG 3079, the GS(2 — 1) line is 
stronger when the beam is offset by +10” in Declination 
(Dec), and separately, when offset by +10” in Right Ascen¬ 
sion (RA). NGG 6946 shows the strongest emission in the 
centre, g etting weaker out towards offsets of ± 20 ” in all di¬ 
rections. Mauersberger et al. (1989) discovered GS in NGG 
6946 and found GS(3 — 2) in offset position (0, 0) and (10, - 
10) from coordinates 0:2000 = 20:34:51 and ^2000 = 60:09:25. 
This is a few arc seconds from our central position. Our re¬ 
sults appear to agree with this but ours also show strong 


^ http://www.iram.fr/IRAMFR/GILDAS 


GS emission from other regions, most notably from the off¬ 
set of +10” in RA. We fit siMle Gaussian profiles to all 
detections of GS (Figure [^. We also detect methanol 
(GH 3 OH) and formaldehjMe (H 2 GO) in both galaxies. The 
isotoplogue of GS, G^^S was marginally detected, also in 
both galaxies. Tables and show the Gaussian fit pa¬ 
rameters for these detections. 

Since we have no prior knowledge as to the location 
of GS emission, we consider each position as separate and 
assume the source fills the beam. The signal-to-noise (S/N) 
of GS in some offset positions is less than 3, and therefore 
they are considered as tentative detections. Still, we used 
them in our analysis, as the fitted line widths and velocity 
positions agree with the values expected. NGG 6946 and 
NGG 3 079 angular sizes are « 12’ x 8 ’ and Ri 4.5’ x 1’ in 
the IR ( Skrutskie et al.|2006 l so observations focus on the 
area surrounding the centre. It is interesting to note that 
in offset positions of just ± 20 ” in any direction in both 
galaxies we fail to see (or only marginally detect) GS. Ball 
et al. (19851 found a bar structure, running approximately 
north-south in NGG 6946 of size 11” by 55”. More material 
in this bar may mean that there is more dense gas, so it 
might be expected that some GS is detected further away 
from the centre, in at least one direction. 


4.2. Line analysis 

Under the assumption that emission from different 
molecules is coming from the same regions, we calculate ve¬ 
locity integrated line ratios for some species. We note that 
this assumption may not hold completely true, but include 
the values as they may be instructive. We calculate the ratio 
of GS(3 — 2) / GS(2 — 1 ) where possible, accounting for the 
difference in beam sizes between the two detections. Since 
we only record detections of lines other than GS(2 — 1) in 
our central pointing in NGG 3079, we only find ratios for 
the central position in this galaxy. These are listed in Ta¬ 
ble 1^ together with ratios in the central position of NGG 
6946. A map of GS(3 — 2) / GS(2 — 1) ratios in NGG 6946 
is shown in Table In NGG 3079 we see much stronger 
emission from GS(3 — 2) when compared to all other lines. 
In the central position, the GS(3 — 2) / GS(2 — 1) ratio 
is approximately double the equivalent in NGG 6946. The 
same is seen when comparing GS(3 — 2) to other lines. In 
NGG 6946 the GS(3 — 2) / GS(2 — 1) ratio is highest at 
offset +10” in RA (approximately 270 pc). Other than this, 
the general trend is for these values to drop as the beam is 
pointed away from the galactic centre, although the value 
is always > 1 . 

Although we have treated lines as being a single peak, 
when the data is less smoothed, detected GS lines in both 
galaxies show a double peak profile (Figure § & [^, al¬ 
though this is not very clear in the tentative detections. 
We consider a possible cause of this to be a rotating torus 
of dense gas around the nucleus, leading to red- and blue- 
shifted comp onents. Profiles like these have been observed 
in NGG 253 (Martin et al. 20051 leading the authors to a 
similar conclusion. However, we note that this might be due 
to other factors, such as self absorption. The lower veloc¬ 
ity (red-shifted) component in both galaxies is narrower, 
possibly meaning that the gas moving away has a lower 
turbulent velocity than the gas coming towards. We show 
the double-horn Gaussian fits for the centre of each galaxy 
and compare with the isotopolgue G^^S in NGG 6946 (Fig- 
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Table 1: Galaxy properties 


Galaxy 
NGG 6946 
NGG 3079 


RA (J2000) Dec (J2000) Velocity (km s Glassification SFR (Mq yr 
20:34:51.88 60:09:14.9 48°^ Scd ^ 

10:01:57.79 55.40:47.0 1124“ SB(s)c 6“ 


NOTES, a) 


Crook et al. 


(20071, b) Leroy et al. i 2008 I, c) 


Mineo et al. 


(20141 





Frequency (MHz) 


Fig. 2: Potential detections in centre of NGG 6946 


ure Detections of isotopologues were marginal (<2 a) 
in NGG 3079 so only the double-horn GS profile is fitted 
(Figure [^. These data have been included with the single 
peak data in Table and For both single and double 
peak profiles, a good agreement can be found between the 
Gaussian fits and the observations. Where possible we fit 
to the double peak but we also calculate for a single fit. 
If a double fit is not possible (as with many of the offset 
points), we only fit a single Gaussian. For double fits we 
have designated the component at lower velocity 1 and the 
higher velocity component 2 (e.g. NGG 6946-1 and NGG 
6946-2). 


4.3. LTE analysis 


We used two methods to perform an LTE analysis of our ob- 
servations. Firstly, we used the rotational diagram method 
( Goldsmith fc Langer]1999 1. With a critical density for low- 
J transitions at low temperature of « 10® cm“® methanol is 
also a high density gas tracer, therefore we use this method 
to constrain a temperature with which to begin an LTE 
analysis using GS {ucrit ~ 1 x 10® for low-J transitions at 
low temperature). We do this only for the centre of NGG 
6946, as lines here are strong enough to allow the use of 
this technique. This method assumes the gas is in LTE and 
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is optically thin. If this is the case then the data are well 
fitted by a linear regression. Under LTE conditions, in the 
rotational diagram method, we are assuming that the ro¬ 
tational temperature is equal to the excitation of all the 
observed transitions. Since the gas may not be thermalised 
the kinetic te mperature may be greater th an the rotational 
temperature (Goldsmith & Langer| 19991. A value for the 
column density can also be obtained using this method. We 
detected two bands of methanol in our spectra (J = 2^ - 
1/c, 3fc - 2fc). Martin et al. (20061 used a method to sepa¬ 
rate individual transitions in these bands using the ratios 
of Einstein coefficients A'^i and A" , upper energy levels 
and E” and degeneracies and g”. This is shown in the 
equation 


W' 

W” 


,,//2 \l _/ 

^ ^ul9u (E':,-E'„)/kT°„, 

7/2 4 " n” ^ 

^ ^ul9u 


( 1 ) 


where W is the integrated flux = J v is the fre¬ 

quency and is the estimate for the rotational tempera¬ 
ture. The method works on the basis that in LTE we should 


see a linear regression in the rotation diagram, the correct 
temperature is the one that produces the best fit. We there¬ 
fore iterate over many temperatures using both methanol 
bands until we find a best fit. The methanol rotation dia¬ 
gram for the centre of NGC 6946 (Figure shows a fit to 
a temperature of 14 K. 

In order to obtain an estimation of the GS column den¬ 
sity, we use only our GS(2 — 1) detections, as we do not 
detect GS(3 — 2) in as many locations as GS(2 — 1). This 
allows us to calculate values in locations where we have 
only one detection. In the same way as with the rotational 
diagram method, we assume LTE conditions and optically 
thin emission. However, it does not constrain a tempera¬ 
ture. A rotational temperature must therefore be assumed. 
Since we obtained a lower limit to the kinetic temperature 
of 14K using methanol, column densities are calculated for 
15 K, 50 K and 300 K. These calculations have been car¬ 
ried out using both single- and (where possible) double-fits. 
They are listed in Table These results are useful as a 
guide to our chemical modelling as they represent a lower 
limit to the column density of GS due to the lack of infor- 
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Table 2: NGC 6946 observational parameters and Gaussian fits. For a two component fit, NGG 6946-1 corresponds to 
the component at lower velocity and NGG 6946-2 the higher velocity component. 


Source 

Line 

V 

Offset 

Offset 

/(T^idjz) 

FWHM 

Tpeafc 

rms 

''^peak 



(GHz) 

RA (") 

Dec (") 

(K km s 

{km s“^) 

(mK) 

(mK) 

(km s“^) 

NGC 6946 

CS(2 - 1) 

97.980 

0 

0 

3.9 ± 0.1 

140.1 ± 4.2 

25.9 

1.6 

60.4 ± 1.7 




0 

-10 

2.2 ± 0.1 

136.6 ± 8.6 

15.0 

2.1 

69.0 ± 4.1 




0 

4-10 

1.9 ± 0.1 

146.8 ± 10.9 

12.2 

2.1 

56.1 ± 4.7 




-10 

0 

2.9 ± 0.1 

139.2 ± 7.3 

19.8 

2.5 

47.0 ± 3.6 




-10 

0 

2.1 ± 0.1 

133.8 ± 7.0 

14.8 

2.0 

70.4 ± 3.5 




-20 

0 

1.6 ± 0.1 

157.6 ± 13.7 

9.3 

1.7 

70.3 ± 6.4 





0 

1.0 ± 0.1 

150.1 ± 18.1 

6.3 

2.2 

33.9 ± 10.3 




-FlO 

4-10 

1.4 ± 0.1 

140.0 ± 15.7 

9.1 

2.3 

41.6 ± 6.2 




+ 10 

-10 

- 

- 

- 

- 

- 




0 

-20 

1.1 ± 0.1 

123.5 ± 15.3 

8.6 

1.9 

71.8 ± 6.6 




0 

+20 

0.9 ± 0.1 

126.0 ± 22.4 

6.4 

2.0 

35.0 ± 10.1 


CS(3 - 2) 

146.969 

0 

0 

3.8 ± 0.1 

135.2 ± 3.7 

26.3 

2.1 

63.3 ± 1.8 




0 

-10 

1.7 ± 0.2 

124.2 ± 16.5 

12.7 

3.3 

74.0 ± 6.7 




0 

4-10 

1.4 ± 0.2 

128.5 ± 15.8 

10.3 

3.9 

54.1 ± 7.7 




-FlO 

0 

3.1 ± 0.2 

114.6 ± 8.8 

25.5 

4.9 

48.1 ± 4.1 




-10 

0 

1.2 ± 0.2 

160.8 ± 7.2 

7.3 

3.4 

79.7 ± 11.4 




-20 

0 

0.7 ± 0.1 

124.3 ± 20.1 

5.3 

2.8 

38.4 ± 11.3 




+20 

0 

- 

- 

- 

- 

- 




+ 10 

4-10 

0.7 ± 0.1 

98.1 ± 19.8 

6.9 

3.5 

35.4 ± 11.0 




-10 

-10 

- 

- 

- 

- 

- 




0 

-20 

- 

- 

- 

- 

- 




0 

+20 

0.4 ± 0.1 

93.6 ± 29.8 

3.9 

3.0 

32.9 ± 15.1 


C'^'‘S(2 - 1) 

96.413 

0 

0 

0.4 ± 0.1 

122.5 ± 25.9 

3.3 

1.4 

65.9 ± 12.7 


C'^'‘S(3 - 2) 

144.617 

0 

0 

0.3 ± 0.1 

68.3 ± 16.2 

3.6 

2.1 

64.7 ± 9.8 


H2CO(2(0,2)-1(0,1)) 

145.603 

0 

0 

1.3 ± 0.1 

161.2 ± 16.4 

7.6 

2.0 

73.0 ± 6.2 


CH 30 H( 2 fc - Ifc) 

96.745 

0 

0 

1.5 ± 0.1 

97.8 ± 6.2 

14.5 

1.6 

71.5 ± 2.5 


CH30H(3fc - 2fc) 

145.125 

0 

0 

2.4 ± 0.1 

106.0 ± 5.0 

21.4 

2.1 

105.3 ± 1.8 

NGC 6946-1 

CS(2 - 1) 

97.980 

0 

0 

1.25 ± 0.2 

65.1 ± 7.2 

18.0 

2.1 

45.1 ± 3.3 



0 

-10 

0.6 ± 0.1 

49.9 ± 10.0 

11.3 

2.1 

-0.8 ± 5.1 




0 

4-10 

- 

- 

- 

- 

- 




-10 

0 

1.23 ± 0.2 

61.1 ± 9.6 

19.1 

2.5 

16.5 ± 5.3 




-10 

0 

0.68 ± 0.2 

59.5 ± 9.4 

10.7 

2.0 

12.8 ± 5.2 




-20 

0 

- 

- 

- 

- 

- 




+20 

0 

0.4 ± 0.1 

42.7 ± 9.5 

9.7 

2.2 

-5.7 ± 5.0 




-bio 

4-10 

- 

- 

- 

- 

- 




+ 10 

-10 

- 

- 

- 

- 

- 




0 

-20 

- 

- 

- 

- 

- 




0 

+20 

- 

- 

- 

- 

- 


CS(3 - 2) 

146.969 

0 

0 

1.2 ± 0.1 

48.1 ± 3.9 

23.0 

2.1 

8.0 ± 1.8 




0 

-10 

- 

- 

- 

- 

- 




0 

4-10 

- 

- 

- 

- 

- 




-FlO 

0 

- 

- 

- 

- 

- 




-10 

0 

- 

- 

- 

- 

- 




-20 

0 

0.4 ± 0.1 

46.6 ± 8.5 

7.9 

2.8 

1.4 ± 5.1 




+20 

0 

- 

- 

- 

- 

- 




+ 10 

4-10 

- 

- 

- 

- 

- 




+ 10 

-10 

- 

- 

- 

- 

- 




0 

-20 

- 

- 

- 

- 

- 




0 

+20 

- 

- 

- 

- 

- 

NGC 6946-2 

CS(2 - 1) 

97.980 

0 

0 

2.5 ± 0.2 

86.7 ± 7.8 

26.6 

2.1 

101.2 ± 4.1 




0 

-10 

1.5 ± 0.1 

80.4 ± 8.2 

17.8 

2.1 

88.0 ± 4.0 




0 

4-10 

- 

- 

- 

- 

- 




+ 10 

0 

1.6 ± 0.3 

73.3 ± 13.5 

19.9 

2.5 

94.9 ± 3.2 




-10 

0 

1.4 ± 0.1 

77.8 ± 9.7 

16.8 

2.0 

96.2 ± 3.5 




-20 

0 

- 

- 

- 

- 

- 




+20 

0 

0.5 ± 0.1 

83.0 ± 19.3 

6.0 

2.2 

85.7 ± 9.9 




-HlO 

4-10 

- 

- 

- 

- 

- 




+ 10 

-10 

- 

- 

- 

- 

- 




0 

-20 

- 

- 

- 

- 

- 




0 

+20 

- 

- 

- 

- 

- 


CS(3 - 2) 

146.969 

0 

0 

2.4 ± 0.1 

79.2 ± 5.1 

28.6 

2.1 

87.1 ± 1.9 




0 

-10 

- 

- 

- 

- 

- 




0 

4-10 

- 

- 

- 

- 

- 




+ 10 

0 

- 

- 

- 

- 

- 




-10 

0 

- 

- 

- 

- 

- 




-20 

0 

0.3 ± 0.1 

54.6 ± 21.4 

5.6 

2.8 

83.4 ± 7.4 




+20 

0 

- 

- 

- 

- 

- 




+ 10 

4-10 

- 

- 

- 

- 

- 




-rlO 

-10 

- 

- 

- 

- 

- 




0 

-20 

- 

- 

- 

- 

- 




0 

+20 

- 

- 

- 

- 

- 


mation about the optical depth. We can therefore discard 
models which show a column density of GS of less than these 
values at a particular temperature. These results show that 


reasonable values of column density (N(GS) cm 

are seen up to temperatures of 300 K. 


Article number, page 6 of 


14 









G. Kelly et al.: Mapping CS in starburst galaxies: Disentangling and characterising dense gas 


Table 3: NGC 3079 observational parameters and Gaussian fits. For a two component fit, NGG 3079-1 corresponds to 
the component at lower velocity and NGG 3079-2 the higher velocity component. 


Source Line v Offset Offset f FWHM "Tp^ak rins ^peak 




(GHz) 

RA (") 

Dec (") 

(K km s 

(km s ^) 

(mK) 

(mK) 

(km s 

NGC 3079 

CS(2 - 1) 

97.980 

0 

0 

2.2 ± 0.2 

429.5 ± 34.6 

4.7 

1.4 

1153 ± 15 




0 

-10 

- 

- 

- 

- 

- 




0 

+ 10 

2.7 ± 0.4 

464 ± 72 

5.4 

2.5 

1213 ± 26 




-fio 

0 

2.4 ± 0.3 

472 ± 67 

4.8 

2.6 

1149 ± 26 




-10 

0 

- 

- 

- 

- 

- 




-10 

+ 10 

- 

- 

- 

- 

- 




-20 

+20 

- 

- 

- 

- 

- 


CS(3 - 2) 

146.969 

0 

0 

4.3 ± 0.3 

484.9 ± 40.8 

8.2 

2.5 

1140 ± 14 


H2CO(2(0,2)-1(0.1)) 

145.603 

0 

0 

1.9 ± 0.2 

444 ± 59 

4.0 

2.7 

1190 ± 28 


CH 30 H( 2 fc - Ik) 

96.745 

0 

0 

0.9 ± 0.1 

362 ± 61 

2.3 

1.1 

1240 ± 31 


CH30H(3fc - 2k) 

145.125 

0 

0 

2.0 ± 0.3 

394 ± 64 

4.7 

3.2 

1158 ± 27 

NGC 3079-1 

CS(2 - 1) 

97.980 

0 

0 

0.9 ± 0.2 

208.3 ± 42.1 

4.3 

1.5 

1010 ± 21 


CS(3 - 2) 

146.969 

0 

0 

2.0 ± 0.4 

244.4 ± 41.1 

7.8 

2.5 

994 ± 22 

NGC 3079-2 

CS(2 - 1) 

97.980 

0 

0 

1.2 ± 0.2 

222.0 ± 35.4 

4.8 

1.5 

1270 ± 20 


CS(3 - 2) 

146.969 

0 

0 

2.1 ± 0.4 

230.3 ± 37.6 

8.5 

2.5 

1279 ± 19 


Table 4: Other line ratios for central pointing 
Galaxy 


Lines 


Ratio 


NGG 6946 GS(3 - 2) / GS(2 - 1) 3.2 ± 0.2 

GH30H(2fc - Ifc) / GS(2-1) 0.4 ± 0.1 

GH30H(3fc-2fc) / GS(3-2) 0.6 ± 0.1 

GH30H(2fc-lfc) / GS(3-2) 0.1 ± 0.006 

GH30H(3fc - 2fc) / GS(2 - 1) 2.0 ± 0.1 

NGG 3079 GS(3 - 2) / CS(2 - 1) 6.7 ± 0.3 

GH30H(2fc - Ifc) / GS(2-1) 0.4 ± 0.1 

GH30H(3fc-2fc) / GS(3-2) 0.5 ± 0.1 

GH30H(2fc-lfc) / GS(3-2) 0.06 ± 0.003 

GH30H(3fc-2fc)/GS(2-1) 3.0 ± 0.2 


Table 5: GS(3 - 2)/GS(2 - 1) ratios in NGG 6946 



RA offset (”) 

20 10 0 -10 

-20 


20 


1.5 


to 

0 

10 

1.7 

2.4 


0 

3.5 

3.2 1.9 

1.5 

P 

-10 


2.5 



-20 





4.4. Non-LTE analysis 
4.4.1. Chemical modelling 

A time-dependent chemical modelling is necessary in or¬ 
der to take into consideration the history of the gas. Here 
we are looking at a dense core environment. To determine 
the properties of the gas in this type of environment, we 
first complete a detailed chemical modelling, and subse¬ 
quently link the outputs of this with a molecular line ra- 
diative transfer model, as detailed in |Benedettini et 
(20061. Since the majority of GS emission is likely to be 


from dense, star forming regions, we model the formation 
of GS through the collapse of diffuse molecular gas to a 
dense core. The chemical model we use is UGL GHEM 


( |Viti et al.|2004 1. UGL_GHEM is a time and depth depen- 
dent gas-grain chemical model operating over two phases. 
Phase I simulates the collapse of initially atomic gas to a 
dense molecular cloud. During this phase, temperature is 
a constant at 10 K. Freeze out of species onto grain sur¬ 
faces occurs; the model uses the chemical network from the 


UMIST database (Woodall et al. 20071 comprising of 211 
species with 2303 gas and grain surface reactions. Phase 
H simulates the gas once it has reached a dynamical equi¬ 
librium. A burst of star formation is simulated, warming 
the cloud, increasing the temperature up to a maximum 
of 400 K. The species included in the model include all 
extragalactic molecules recorded so far, as well as many 
galactic species. Surface reactions are mainly based around 
hydrogenation, allowing saturated species to form. We run 
a large grid of models varying four parameters. These are, 
the temperature in phase H (T), the UV radiation field (%), 
the cosmic ray ionisation rate (C) and the final gas density 
of the cloud in phase I (n(H 2 )). Each model is ran over a pe¬ 
riod of 10® years. All other parameters are taken as solar, 
scaled to the metallicity of the two galaxies where neces¬ 
sary. We ran a large grid of models (summarised in Table 
[ 7 1, varying parameters over a wide range, given the lack 
OT prior knowledge of the environment we are considering. 
By doing this we can see the effect on the chemistry and 
ultimately the formation of GS in different environments. 


4.4.2. Radiative transfer modelling 

The chemical abundances of GS given by UGL GH EM are 
are input to the radiative transfer model SMMOL (Rawl¬ 


ings & Yates 20011. This is an accelerated A-iteration (ALl) 
code for solving multi-level radiative transfer problems. 
Data on the GS molecule - including energy levels, tran¬ 
sition frequencies, Einstein A coefficients and collisional 


rates - are taken from the LAMDA database (Schoier et al. 


20051. Initially the code calculates the radiation held and 
the level populations. This is completed under the assump¬ 
tions of LTE cond itions and ISM radiation as input contin¬ 
uum ( Black[19M I. The radiation field is then recalculated, a 
check is made for convergence. This repeats until the model 
converges. At each radial point the code generates the ra¬ 
diation field and the level populations. The emergent in¬ 
tensity distributions are then convolved with the telescope 
beam. This means the model directly predicts the line pro- 
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Fig. 4: Single component fits for CS(2 — 1) and CS(3 — 2) in NGC 6946. Label in right corner shows offset in RA (”) and 
Dec (”) from central position. 


file for a given source as observed with a given telescope. 
The code has been successfully benchmarked against other 
radiative transfer models. The ability to predict a line pro¬ 
file allows us to compare with our observations, giving us a 
single model that best fits the chemical history and physi¬ 
cal conditions of the environment observed. Both SMMOL 
and UCL_CHEM are modelling the formation of one dense 
core. In our beam, there must be many more dense cores. 
At most, the beam covers an area ^2 kpc (CS(2 — 1) in 
NGC 3079) and the typical scale of dense cores are 0.05 


pc. It is extremely unlikely CS emission is emanating only 
from one dense core. If it was, we would not see it as the 
signal would be too diluted in the beam. It is more likely 
that the emission is coming from a number of these dense 
core regions. To account for this, we estimate the number 
of cores in our beam and scale up the emergent flux by this 
factor . Our estimate is based on a study by [Lintott et al.| 
(2005|). They estimate that there are 10^ hot dense cores in 
the Milky Way, and calculate a scale factor for other galax¬ 
ies based on their star formation rate. We adopt the same 
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Fig. 5: Single component fits for CS(2 — 1) and CS(3 — 2) in NGC 3079. Label in right corner shows offset in RA (”) and 
Dec (”) from central position. 


methodology, based on the star formation rates of our ob¬ 
served galaxies (see Table and the likely central location 
of the gas. Our targets are both starburst galaxies which 


concentrated over the nuclear bulge of both galaxies (Tsai 

et al. 

2013 

Yamagishi et al. 

2010 

Schinnerer et al. 

2006p. 


ders of magnitude, and by the same values for both galaxies, 
although we acknowledge that the regions observed in NGC 
3079 are larger than in NGC 6946. We use scale factors of 


1 , 10, and 100; leading to our estimates for the number of 
dense cores to be 10^, 10®, and 10®. We note that we take 
the arbitrary end point of the chemical model (10® years) 
to input into SMMOL. While we do not claim this to be the 
age of our cores, we use it as a representative chemical age 
of the molecular gas at large scales. In fact, time dependent 
effects may be important in energetic extra galactic environ¬ 
ments and have been dis cussed elsewhere (Meijerink et al. 
2013 Bayet et al.||2008b l. 
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Table 6 : LTE Column Densities 


Source 

Line 

Offset 

Offset 


N(CS) X 10^^ 




RA (") 

Dec (") 

T = 15 K 

T = 50 K T 

= 300 K 

NGC 6946 

CS(2- 1) 

0 

0 

0.176 

0.42 

2.25 


0 

-10 

0.099 

0.24 

1.27 



0 

+10 

0.087 

0.21 

1.11 



+ 10 

0 

0.133 

0.32 

1.71 



-10 

0 

0.096 

0.23 

1.22 



-20 

0 

0.071 

0.17 

0.91 



+20 

0 

0.046 

0.11 

0.58 



+ 10 

+10 

0.062 

0.15 

0.79 



+ 10 

-10 

- 

- 

- 



0 

-20 

0.051 

0.12 

0.66 



0 

+20 

0.039 

0.09 

0.5 

NGC 6946-1 

CS(2- 1) 

0 

0 

0.057 

0.14 

0.73 


0 

-10 

0.027 

0.07 

0.35 



0 

+10 

- 

- 

- 



+ 10 

0 

0.056 

0.14 

0.72 



-10 

0 

0.031 

0.07 

0.4 



-20 

0 

- 

- 

- 



+20 

0 

0.02 

0.05 

0.26 



+ 10 

+10 

- 

- 

- 



+ 10 

-10 

- 

- 

- 



0 

-20 

- 

- 

- 



0 

+20 

- 

- 

- 

NGC 6946-2 

CS(2- 1) 

0 

0 

0.111 

0.27 

1.43 


0 

-10 

0.069 

0.17 

0.89 



0 

+10 

- 

- 

- 



+ 10 

0 

0.071 

0.17 

0.9 



-10 

0 

0.063 

0.15 

0.81 



-20 

0 

- 

- 

- 



+20 

0 

0.024 

0.06 

0.31 



+ 10 

+10 

- 

- 

- 



+ 10 

-10 

- 

- 

- 



0 

-20 

- 

- 

- 



0 

+20 

- 

- 

- 

NGC 3079 

CS(2- 1) 

0 

0 

0.098 

0.24 

1.26 



0 

-10 

0.031 

0.08 

0.4 



0 

+10 

0.122 

0.29 

1.56 



+ 10 

0 

0.11 

0.26 

1.41 



-10 

0 

0.067 

0.16 

0.86 



-10 

+10 

0.074 

0.18 

0.95 



-10 

+20 

0.045 

0.11 

0.58 


indicates no fit was possible 


Table 7: All permutations of these parameters were ran 
through UCL_CHEM 

T (K) X (Habing) C (1-3 x IQ-^^ s"^) n(H 2 ) (cm-^ 


30 

1 

1 

10 ^ 

50 

10 

10 

10 “^ 

100 

100 

100 

10 ^ 

200 

1000 

1000 

10 ® 

300 

10000 

10000 

10 ^ 

400 

- 

- 

10 ® 


Analysis of the models is completed using a method. 
We use the equation 


where N is the number of lines used, Thiod is the integrated 
flux of the modelled line, and Fobs is the integrated flux of 
the observed line.This form ula has been used in previous 
papers in a similar way (e.g. Benedettini et alT]|2006 Doty 


et al.|2004). We use the difference between the model flux 


and the observed flux instead of the measurement error be¬ 
cause the errors associated with the modelling are greater 
than those in the observations. This also has the advan¬ 
tage of avoiding giving brighter lines a higher weight. It is 
possible to use SMMOL to model line profiles of CH 3 OH 
and H 2 CO. Since neither of these molecules were detected 
outside of the centre of either galaxy, it was not possible in 
this study. 


Article (i) 

^ Fobs(*) 


1 2 




(2) 
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Table 8: Best fit parameters for each location in each galaxy (10® cores) 


Galaxy 

Aa 

A6 

T (K) X (Habing) C ( Co = 1-7 x lO"!" s'^) 

n (cm ®) 


Notes 

NGC 6946 

0 

0 

100 1 

100 

10 ® 

0.20 


NGC 6946 

0 

-10 

100 1 

1 

10 ® 

0.032 


NGC 6946 

0 

+ 10 

100 1 

1 

10 ® 

0.046 


NGC 6946 

+10 

0 

100 1 

1 

10 ® 

0.021 


NGC 6946 

-10 

0 

100 1 

1 

10 ® 

0.13 


NGC 6946 

-20 

0 

100 1 

10 

10 " 

0.55 

poor fit 

NGC 6946 

+20 

0 

300 1 

100 

10 " 

0.06 

(2 - 1) only 

NGC 6946 

+10 

+ 10 

100 1 

10 

10 " 

1.1 

poor fit 

NGC 6946 

+10 

-10 





no fit 

NGC 6946 

0 

-20 

400 1 

1000 

10 " 

0.02 

(2 - 1) only 

NGC 6946 

0 

+20 

100 1 

10 

10 " 

0.80 

poor fit 

NGC 3079 

0 

0 

400 1 

10 

10 ® 

0.088 


NGC 3079 

+10 

0 

400 1 

10 

10 ® 

0.009 

(2 - 1) only 

NGC 3079 

0 

+ 10 

400 1 

10 

10 ® 

0.08 

(2 - 1) only 




Table 9 

: Best fit parameters for each location in 

each galaxy (10® 

cores) 


Galaxy 

Aa 

A,5 

T (K) 

X (Habing) 

C ( Co = 1.7 X 10- 

"" s”") n (cm” 

") X^ 

Notes 

NGC 6946 

0 

0 

400 

1 

10 

10 ® 

0.30 


NGC 6946 

0 

-10 

200 

1 

1 

10 ® 

0.24 


NGC 6946 

0 

+10 

100 

1 

10 

10 ® 

0.25 


NGC 6946 

+ 10 

0 

400 

1 

10 

10 ® 

0.22 


NGC 6946 

-10 

0 

100 

1 

10 

10 ® 

0.24 


NGC 6946 

-20 

0 

100 

1 

10 

10 ® 

0.21 


NGC 6946 

+20 

0 

400 

1 

10 

10 ® 

0.0001 

(2-1) only 

NGC 6946 

+ 10 

+10 

100 

1 

10 

10 ® 

0.31 


NGC 6946 

+ 10 

-10 






no fit 

NGC 6946 

0 

-20 

400 

1 

10 

10 ® 

0.01 

(2-1) only 

NGC 6946 

0 

+20 

400 

1 

100 

10 ® 

0.02 

(2-1) only 

NGC 3079 

0 

0 

400 

1 

1 

Io®“ 

0.37 


NGC 3079 

+ 10 

0 

400 

1 

1 

10 ® 

0.02 

(2-1) only 

NGC 3079 

0 

+10 

400 

1 

1 

10 ® 

0.05 

(2-1) only 


5. Discussion 


We now analyse the results of our models. It is important 
to note that while the observations cover a large area, CS 
emission is predominantly emanating from smaller regions 
of dense, star forming gas. We assume that there are many 
of the s e regi ons in our beam, based on the work by |Levine| 
et al. (20081. Since we complete this study under the as¬ 
sumptions either 10"*^, 10® or 10® dense cores in the beam, 
we discuss each set of models separately. Firstly, under the 
assumption of 10® cores in the beam, the best fit models for 
all positions in NGC 3079 are the same (Table |^. Tem¬ 
perature and density are quite high - 400 K & 10° cm“®, 
respectively - with a cosmic ray ionisation rate of 100 times 
the standard rate (Co). The high temperature is likely to be 
somewhat higher than the average temperature across the 
region we are observing. The high temperature is possibly 
an indication of the presence of shocks, which we discuss 
later. Models for NGC 6946 follow the general pattern of 
decreasing density (10® cm“® —)■ 10® cm“®) moving away 
from the centre of the galaxy (when only including models 
with good fits to both lines). Models for the outer regions of 
our observations generally do not provide good fits, or there 
is not enough data to fit to. Although the lowest x^s are 


all for high density (10^ cm“®) models. We include these 
results for completeness, but note that - given the high 
values - there are large degeneracies in these models. There 
are no models with a density <10® cm“® that provide a 
good fit to the data (i.e. with <!)• 

When considering 10® cores in the beam, the best fit 
models for each position (see Table in NGC 3079 arise 
from a different set of parameters from the previous set 
of models, although again there is no variation across posi¬ 
tion. Again from a high temperature model, the density and 
cosmic ray ionisation rate are lower than for the 10® core 
simulation but do not provide the same quality fit as shown 
in Table [S] In contrast with the results for NGC 6946 in 
Table the 10® cores models show a higher temperature 
and higher density for some best fit models, although as 
with NGC 3079, these do not fit as well as the 10® core 
models. The UV radiation field has no effect on any of the 
best fit models for any number of cores, because of the high 
visual extinction associated with the high density regions 
where CS emission originates. We note that this does not 
imply there is a lack of UV radiation but simply that UV 
is unable to penetrate down to dense cores. We also find 
that we cannot successfully fit the observations using 10^ 
cores for any region of either galaxy because the antenna 
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Fig. 6: CS detections in centre of NGC 6946. dis¬ 

playing a double-horn profile. 



500 1000 1500 2000 

Velocity (km/s) 


30 

29.5 

29 

28.5 

^ 28 

27.5 

27 

26.5 

26 

5 10 15 20 25 30 35 40 45 50 55 

Eu (K) 

Fig. 8: Rotation diagram for methanol in centre of NGC 
6946. The parameters for this diagram are T = 13.6 K and 
N(CH30H) = 4.6 X cm-2. 
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"H, 


'Q 


20011. We ran shock models for the central position of both 
galaxies. These were completed with the same best fit pa¬ 
rameters as in the 10® core models, as these give the low¬ 
est values. During the shocked period, the temperature 
rises from 100 K to 500 K, and then reduces back to 100 K. 
During this period, anything remaining on grain surfaces 
is sublimated. In our radiative transfer code we then vary 
the number of cores as before. We see that the best fits to 
the line intensities of CS(2 — 1) and CS(3 — 2) now come 
from the model assuming that there are 10^ dense cores 
within the beam (see Figures [^& [Io| . This suggests that 
since shocks have a significant effect on CS emission, it is 
necessary to know whether the region being observed is (or 
has been) shocked, in order to fully characterise the envi¬ 
ronment observed. Therefore, observations of shock tracers 
(such as SiO), would be very useful in conjunction with CS 
detections. 

The line profiles calculated by SMMOL are not just de¬ 
pendent on the physical conditions of the environment, but 
also on the chemistry of CS. In the following section we 
analyse the outputs from UCL_CHEM in order to deter¬ 
mine what factors influence CS chemistry. 


Fig. 7: CS detections in centre of NGC 3079. dis¬ 

playing a double-horn profile. 


temperature of the theoretical lines is too low to match ob¬ 
servations. In many cases, the profiles do not match either. 
This possibly provides a lower limit on the number of star 
forming cores we see within our beam. 


5.2. CS formation 


The dominant formation routes to CS are not com pletely 
understood. There was initially some agreement (Drdla 


et al.|19^ Lucas fc Liszt|2002| that most CS formation in 


the gas phase is a product of the dissociative recombination 
reaction: 


HCS+ -f e 


CS-kH, 


(3) 


5.1. Shocks 

All the discussion above assumes there are no shocks. How- 
CS has recently been shown to be present in shocked 


ever, 


regions (e.g. Nicholas et al. 


the high temperature models which successfully fit many of 
the emission regions. Shock temperatures are unlikely to be 
sustained for the timescale of our models but could be very 
significant for CS chemistry even over a shorter timescale. 
UCL_CHEM is capable of simulating shocks (Viti et al. 


What is less clear is the route to the formation of HCS'*'. 


20121. This is also indicated by 


Drdla et al. (19891 found that the following gas-phase re- 
action scheme reproduced CS in abundances in the same 
range as their observations in relatively diffuse regions: 


S+ + 
S+ + 
CS+ 


CH - 
C2- 
FH2 


> CS+ -k H, 
CS+ -kC, 
HCS+ -k H, 


( 4 ) 
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Fig. 9: Shock model fits for NGC 6946 Fig. 10: Shock model fits for NGC 3079 


However, Lucas & Liszt (20021 found that the abun¬ 
dances of S+, G 2 , and GH were an order of magnitude 
lower than observed if the mechanism above was producing 
the observed quantities of GS. Its high efficiency at low 
temperatures mean that reaction is still likely to 
contribute to GS formation in certain environments. A 


general, more well agreed upon formation route (Destree 

et al. 

2009 Nilsson et al. 

to 

0 

0 

0 

Lucas & Liszt||2002p is the 


SO-hC 


CS-hO, 


(5) 


Ghemical modelling ( Bergin fc Langei]|1997 ) has shown 
that this reaction leads to a high abundance of GS quite 
swiftly, until formation of GO locks away much of the 
atomic carbon needed for reaction . At later times, GS 
abundance drops as a result of efficient destruction of GS 
with atomic oxygen: 


CS-hO 


CO-hS, 


( 6 ) 


Routes |5]and |6]are now believed to dominate the forma¬ 
tion and destruction of GS ( Lucas fc Liszt|[2002 1 . Analysis 
of our chemical models tend to broadly a,gree with these 
findings. However, the cosmic ray ionisation rate is an im¬ 
portant parameter in determining the dominating mecha¬ 
nisms. A high cosmic ray ionisation rate (such as the one 


in one of our best fit models for the centre of NGG 6946 - 
C = 100 Co) leads to the destruction of a small percentage 
of available GO, providing enough atomic carbon to allow 
reaction [5] to counteract the destruction from reaction [6l 
However, this is also more efficient as there is more atomic 
oxygen available from GO. At a lower cosmic ray ionisation 
rate, we find that a significant source of GS destruction is 
from reaction Both the formation and destruction pro¬ 
cesses are up to two orders of magnitude slower when C is 
lowered from 100 Co to 1 Co- 


H30+-fGS—^HCS+-bH20, (7) 

During shock treatment the following formation reac¬ 
tion also becomes important: 


C 2 -f S —^ GS -b C. (8) 

Reaction is important during the shock as it both 
becomes more efficient at high temperatures and the abun¬ 
dance of the reactant species increases. This points to there 
not being one overwhelmingly dominant way of forming GS. 
Different reactions are responsible for its formation (and 
destruction) in different conditions. 

6. Conclusions 

GS emission has been mapped across two nearby starburst 
galaxies - NGG 3079 and NGG 6946. LTE calculations were 
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undertaken to determine a lower limit to the kinetic tem¬ 
perature in NGC 6946 of 14K. In order to constrain this 
further, and determine conditions in the dense gas envi¬ 
ronments being traced by CS, a time dependent gas-grain 
chemical model (UCL_CHEM) has been linked to a radia¬ 
tive transfer code (SMMOL) to simulate the line emission. 
Comparing the observations to the models has led to the 
conclusion that, given the high temperatures (T = 400 K) 
often needed for successful fits, models featuring the pres¬ 
ence of shocks best reproduce some of our observations. We 
find that in order to match the observations, emission must 
be coming from between lO"* and 10® dense cores of ~ 0.05 
pc. The routes to the formation of CS are important fac¬ 
tors in successfully modelling regions of dense gas. With 
ever more complex chemical webs, predicted abundances 
of molecules become highly sensitive to set conditions. We 
find that shocks provide conditions that are effective in effi¬ 
ciently forming CS in necessary abundances due in part to 
increased efficiency of Reaction A high cosmic ray ioni¬ 
sation rate (~100 Co) could also have a significant influence 
on CS chemistry by increasing the rate of Reaction ^ In 
conclusion, CS observations, combined with chemical and 
radiative transfer modelling are capable of disentangling 
dense star forming regions from surrounding molecular gas. 
With better spatial resolution it should be possible to quan¬ 
tify further the chemical differentiation across the central 
regions of starburst galaxies. Future multi-line observations 
of CS with ALMA would provide the tools to do this. 
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